arXiv:1507.06991vl [cond-mat.str-el] 24Jul2015 


Fate of the Kondo Effect and Impurity Quantum Phase Transitions Through 

the Lens of Fidelity Susceptibility 


Lei Wang 1 , Hiroshi Shinaoka 1,2 and Matthias Troyer 1 
1 Theoretische Physik, ETH Zurich, 8093 Zurich, Switzerland and 
2 Department of Physics, University of Fribourg, 1700 Fribourg, Switzerland 

The Kondo effect is an ubiquitous phenomenon appearing at low temperature in quantum confined 
systems coupled to a continuous bath. Efforts in understanding and controlling it have triggered 
important developments across several disciplines of condensed matter physics. A recurring pattern 
in these studies is that the suppression of the Kondo effect often results in intriguing physical 
phenomena such as impurity quantum phase transitions or non-Fermi-liquid behavior. We show 
that the fidelity susceptibility is a sensitive indicator for such phenomena because it quantifies the 
sensitivity of the system’s state with respect to its coupling to the bath. We demonstrate the power 
of fidelity susceptibility approach by using it to identify the crossover and quantum phase transitions 
in the one and two impurity Anderson models. 
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The Kondo effect [1] was first observed in 1934 [2] 
as a low temperature resistance minimum in gold and 
was explained by Kondo in 1964 by taking into ac¬ 
count the scattering of conduction electrons and mag¬ 
netic impurities [3]. However, Kondo’s perturbative cal¬ 
culation exhibits an unphysical divergence of the resis¬ 
tance at zero temperature. Resolving the Kondo prob¬ 
lem has ultimately led to significant theoretical pro¬ 
gresses, including the formulation of the scaling laws [4], 
the development of the numerical renormalization group 
(NRG) method [5] and the application of phenomenologi¬ 
cal Fermi-liquid theory [6], Bethe ansatz [7, 8] and bound¬ 
ary conformal field theory [9] to the quantum impurity 
problems. Experimental interests has increased in the 
late 1990s due to breakthroughs in fabricating artificial 
nano-devices and creating tunable Kondo effects [10-15]. 
Other relevant experimental systems include the dissipa¬ 
tive two-state systems [16] and the heavy-fermion com¬ 
pounds [17, 18]. Moreover, through the dynamical-mean- 
field-theory (DMFT) framework [19, 20] a connection be¬ 
tween quantum impurity problems and correlated lattice 
models has been established [21, 22]. 

A general description of the quantum impurity prob¬ 
lems can be written as 

H( A) = -Himpurity + -Hbath +Ai7l, (1) 

S -N,-' 

Ho 


where Hq describes the quantum impurity together with 
a continuous bath, and the last term describes the cou¬ 
pling between them. We treat A as a parameter and aim 
to characterize the state of the quantum impurity as a 
function of its coupling to the bath. The Kondo effect 
originates from the bath’s tendency to screen the local 
moment formed on the quantum impurity. Renormaliza¬ 
tion group analysis shows that in the Kondo region the 
coupling strength flows to infinity at low energy [4, 5] , im¬ 
plying that the local moment will eventually get screened 


at low enough temperature even with an arbitrarily weak 
bare impurity-bath coupling strength. 

There are, however, various physical processes that can 
compete with the Kondo effect. In the presence of such 
competitions the system may undergo an impurity quan¬ 
tum phase transition where a competing state (local mo¬ 
ment, charge order etc.) takes over as the bath-impurity 
coupling A decreases. Suppression of the Kondo screening 
often leads to non-Fermi liquid behavior [23, 24]. How¬ 
ever, different from the quantum phase transition in bulk 
systems [25], at such impurity quantum critical point only 
a non-extensive term in the free energy becomes singu¬ 
lar. It is not always straightforward to find local probes 
to identify the impurity phase transitions. The question 
arises how to diagnose and characterize such impurity 
quantum phase transitions in a general setting [26] . 

In this Letter we argue that the fidelity susceptibil¬ 
ity [27, 28] provides a general and direct probe for an 
impurity quantum phase transitions. With its origin in 
quantum information and in the differential geometry 
perspective of quantum states, the fidelity susceptibility 
does not depend on details of the physical systems. The 
quantum fidelity F(Ai,A 2 ) is defined as the overlap of 
two ground state wave-functions (or density matrices in 
the nonzero temperature case [29]) for coupling strengths 
Ai and A 2 . The fidelity susceptibility [27, 28] 


Xf( A) 


d 2 lnf(A,A + e) 
de 2 


e=0 


( 2 ) 


typically exhibits a maximum at the phase boundary be¬ 
cause the system’s state changes drastically around the 
quantum critical point. Since the fidelity susceptibility 
also fulfills the scaling laws [28, 30], it is an effective tool 
to detect and characterize various quantum phase tran¬ 
sitions, see Ref. [31] for a review. 

Recently, some of us developed an efficient approach 
for calculating the fidelity susceptibility of quantum 
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many-body systems [32] using modern quantum Monte 
Carlo (QMC) methods [33-43] . Specializing this to quan¬ 
tum impurity models one can perform an expansion for 
the bath-impurity coupling [44] 


z = J2 xk dri ■ ■ ■ 


dr k x 


k =0 


Tr 
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Depending on details of the impurity Hamiltonian, vari¬ 
ous QMC algorithms can be used to sample Eq. (3). For 
example, the continuous-time hybridization expansion al¬ 
gorithm (CT-HYB) [39] solves (multi-orbital) Anderson 
impurity models, while the CT-J algorithm [45] is suit¬ 
able for the Kondo model and its multi-orbital general¬ 
ization such as the Coqblin-Schrieffer model [46]. 

Equation (3) also provides a conceptual framework 
to understand an impurity quantum phase transition 
through a quantum-classical mapping. The expansion 
can be formally interpreted as a grand canonical partition 
function of classical particles residing on a ring of length 
/3. These particles represent the bath-impurity coupling 
events provided by the Hi terms and their number is 
controlled by the coupling strength A. Since Eq. (3) has 
the form of a fugacity expansion, an impurity quantum 
phase transition driven by A will manifest itself as a con¬ 
densation phase transition of classical particles [47, 48] . 

A concrete example of this general reasoning is pro¬ 
vided by the Anderson-Yuval solution of the anisotropic 
Kondo model [49-52] with Hamiltonian 

-HKondo = ^ ' e kCj ia .Ck(T T J z ^ ^ S ^cj^Ckq 

k,<T k,k' 

+A ('S ,+ 4;Ck't + h.c) . (4) 

k.k' 

The last term, which plays the role of XHi in Eq. (1) 
describes the coupling of the local impurity spin to free 
electrons in the bath through spin-flips. An expansion 
in the form of Eq. (3) and integration out of the free 
fermions lead to Anderson and Yuval’s mapping of the 
Kondo model to a one-dimensional classical Coulomb 
gas [49]. As illustrated in Fig. 1(a), the spin flips can be 
interpreted as alternating positive and negative charges 
interacting through a logarithmic Coulomb potential [49— 
52]. The Coulomb gas picture provides an intuitive un¬ 
derstanding of the Kondo effect and the renormalization 
group flow [50] . For a ferromagnetic coupling J z < 0, the 
Coulomb gas exhibits a phase transition as the fugacity 
A changes. For small A < | J. | these Coulomb charges are 
dilute and all associated in pairs, corresponding to the 
ferromagnetic Kondo state where the quantum impurity 
is spin polarized, while A > \J Z \ corresponds to the an¬ 
tiferromagnetic Kondo state where the spin-flips are so 
frequent that the impurity shows no net magnetization, 
i.e. it is Kondo screened [53]. 
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Figure 1. (color online), (a) The Anderson-Yuval mapping 
of the Kondo model (Eq. (4)) to a one-dimensional classical 
Coulomb gas. The magnetization of the quantum impurity 
flips in the imaginary time due to the coupling to the bath. 
The spin flips can be interpreted as positive (red circle) and 
negative (blue circle) charges distributed on a periodic ring, 
(b) The continuous-time hybridization-expansion QMC algo¬ 
rithm, in the same spirit, maps the single impurity Anderson 
model (Eq. (6)) to a classical statistical problem. The thick 
segments indicate the occupation of the spin up and down 
impurity levels. The endpoints of each segment represent the 
hybridization events where the electron hops in/out of the 
impurity, which are treated as classical objects in the QMC 
sampling. In the both cases (a) and (b), the fidelity sus¬ 
ceptibility is calculated as the covariance of fcn and which 
count the number of bath-impurity coupling events in the two 
equal-bipartitions of the imaginary-time axis. 


In the framework of Eq. (3), the fidelity susceptibility 
(2) can be readily calculated using a covariance estima¬ 
tor [32] 


XF 


( k L k R ) - (k L ) ( k R ) 
2A 2 


( 5 ) 


where k l and k r are the numbers of Hi operators in 
the two bipartitions of the imaginary-time axis. In the 
case of the Kondo model Eq. (4), they correspond to 
the number of charges on either side of the bipartition, 
shown in the bottom of Fig. 1(a). It is clear from Ander¬ 
son and Yuval’s classical Coulomb gas picture that the 
fidelity susceptibility estimator (5) captures the critical 
fluctuation upon a condensation phase transition there¬ 
fore is able to signify the impurity phase transitions of 
the anisotropic Kondo model. 

For general quantum impurity models the estimator 
(5) always quantifies the sensibility of the system’s state 
with respect to the bath-impurity coupling and can there¬ 
fore be used to diagnose impurity quantum phase tran¬ 
sitions. It is a generic probe of phase transition irre¬ 
spective of physical details of the system. The singular¬ 
ity of the fidelity susceptibility upon a phase transition 
is also stronger than the second order derivative of the 
free energy (related to variance of the total expansion 
order) [30, 32], Moreover, the fidelity susceptibility can 
also be used to inspect the crossover physics, even though 
there is no sharp phase boundary. 
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As illustration we consider first the single impurity An¬ 
derson model [54] 

-ffsiAM = ^2 £ kcj ia c k<7 + e d ^2 + u 

k.er <7 

+ X Yl ( 6 k<T^ + h - C ) > ( 6 ) 

k ,a 

where n a = d\,d a is the impurity occupation number and 
the second line describes the hybridization of the impu¬ 
rity and the noninteracting bath with strength A. We 
consider a noninteracting bath with semicircle density- 
of-states p(e) = ]C k <5(e “ £ k) = ^dV 1 - ( e / D ) 2 with 
D = 2 and choose ed = —U/2 such that model is at 
the particle-hole symmetric point. As we tune the on¬ 
site interaction U and the hybridization strength A there 
is a crossover from a local moment regime, where the 
spin of the singly occupied impurity is free to flip to the 
Kondo region, where the local moment is screened by the 
bath [55]. 

We use the CT-HYB algorithm [39] for our simulations 
and illustrate one Monte Carlo configuration in Fig. 1(b). 
Each dashed line indicates a hybridization event, where 
the electron hops on or off the impurity site, thus chang¬ 
ing the occupation (indicated by the thickness of the seg¬ 
ments). The fidelity susceptibility Eq. (5) is easily mea¬ 
sured by counting the number of hybridization events in 
a bipartition of the imaginary time axis, shown in the 
bottom of Fig. 1(b). 

Figure 2(a) shows the fidelity susceptibility \f( A) in 
the U — A plane with fixed inverse temperature /3 = 100, 
where the peak indicates the crossover from local mo¬ 
ment to Kondo region. This figure is a two-dimensional 
slice of the phase diagram of the Anderson impurity 
model sketched in the seminal NRG work of Ref. [55] 
(Fig. 12). The red solid line shows the contour deter¬ 
mined by the Kondo temperature 1//3 = Tk(U, A) = 
A \/Ue~ 7TU ^ 8X ) [56]. This boundary agrees with the 
maxima of the fidelity susceptibility, showing that it in¬ 
deed correctly captures the crossover physics. The peak 
of fidelity susceptibility is higher at small A region, which 
is a manifestation of the Anderson orthogonality catas¬ 
trophe [57]: even a weak coupling to the quantum impu¬ 
rity drastically change the state of the system. 

To further confirm the relevance of the peak of the 
fidelity susceptibility we calculate the local spin suscep¬ 
tibility 


Xs = f dr{S z {r)S z { 0)), (7) 

Jo 

where S z = (n-j-—rij,)/2 is the magnetization on the impu¬ 
rity. Figure 2(b) shows that 4Ty s , which corresponds to 
the effective moment on the impurity, changes from one 
in the local moment region to zero in the Kondo region. 
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Figure 2. (color online). Crossover from the local mo¬ 
ment to the Kondo region in the single impurity Anderson 
model (Eq. (6)) revealed by (a) fidelity susceptibility (b) spin 
susceptibility 4:T\s defined in Eq. (7). The red solid line 
shows the contour determined from the Kondo temperature 
1//3 = Tk{U, A) [55, 56]. 


The crossover region agrees with the peak determined 
from the fidelity susceptibility in Fig. 2(a). 

As a second example we consider the two-impurity An¬ 
derson model [58] 

7?TIAM — ^ Cd ^ ) Uac 7 Y U ^ 

k,a,cr ex.,a a 

+ JSi ■ $2 + A ^2 (ckaadota + h.C.^j , (8) 

k,a,cr 

where a = {1,2} labels two impurity sites with occu¬ 
pation number h aa = d^daa- The impurities have the 
same local interaction U and onsite energy ed = —U/2. 
Each impurity is coupled to its own bath with the hy¬ 
bridization strength A. The last term represents the 
Ruderman-Kittel-Kasuya-Yosida (RKKY) [59-61] inter¬ 
action between magnetic impurities in a metal. In the ab¬ 
sence of this term, each impurity is Kondo screened by its 
own bath for the choice of A = 1 and j3 = 100. However, 
the antiferromagnetic RKKY coupling J > 0 favors a sin¬ 
glet formed between the two impurity spins, which com¬ 
petes with the Kondo screening and causes an impurity 
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Figure 3. (a) The fidelity susceptibility (b) equal-time spin- 
spin correlation (c) density-of-states at the Fermi level of the 
two-impurity Anderson model Eq. (8) as a function of the 
inter-impurity RKKY coupling strength. The dashed vertical 
line indicates the critical point J c = 0.37 where the equal-time 
spin-spin correlation 3 (iSfSf) = —0.25 [63, 67]. 

quantum phase transition [62, 63]. Detailed studies of 
the two impurity Anderson (and Kondo) model have pro¬ 
vided insights into various aspects of the Kondo [64, 65] 
and heavy fermion physics [66]. 

The simulation of this model goes beyond the segment 
picture illustrated in Fig. 1(b), and we thus adopt an al¬ 
gorithm [44, 68] suitable for general interactions. The fi¬ 
delity susceptibility is still calculated in the same way, by 
simply counting the number of hybridization events. As 
shown in Fig. 3(a), it exhibits an increasingly sharp peak 
as the inverse temperature /3 increases. The peak location 
shifts towards the vertical dashed line, where the equal¬ 
time spin-spin correlation (Si ■ S 2 ) = 3(SfiS|) = —0.25 
in Fig. 3(b). According to previous NRG studies [63, 67] 
the quantum critical point is right at the dashed line. 
Obviously, the fidelity susceptibility offers a better indi¬ 
cation of the phase transition compared to the spin-spin 
correlations because the later quantity is featureless at 
the critical point and has much weaker temperature de¬ 
pendence. Figure 3(c) shows the density of states at the 
Fermi level, which decreases as the spin singlet state takes 
over the Kondo state in the large J limit [69]. 

We next perform a scaling analysis of the fidelity sus¬ 
ceptibility close to the quantum critical point [30]. Since 
an infinite bath was assumed, the only finite dimen¬ 
sion is the inverse temperature. Figure 4 shows the 
scaled fidelity susceptibility Xf/P versus (J — J c )j3 1 ^ 2 
with J c = 0.37, which results in a good data collapse. Al¬ 
though the scaling form is chosen empirically according 



Figure 4. (color online). Data collapse of the scaled fidelity 
susceptibility with J c = 0.37. The data are the same as the 
one in Fig. 3(a). 

to the one for lattice systems [30], the observed scaling 
exponents agree with the considerations of Ref. 70. The 
observed data collapse suggests that the fidelity suscep¬ 
tibility not only captures the impurity quantum critical 
point but also the values of the critical exponents, which 
are an indicator for the universality class of a quantum 
phase transitions. 

Our paper shows that the fidelity susceptibility is a 
versatile tool to probe and inspect phase transition and 
crossover physics in quantum impurity models. It is read¬ 
ily acccessible in QMC simulations and it serves as a gen¬ 
eral purpose indicator for the breakdown of the Kondo 
effect. Conceptually our work exploits the intrinsic quan¬ 
tum to classical mapping of the quantum impurity mod¬ 
els in the context of modern QMC approaches. 

Recent experimental and theoretical studies explore an 
even richer variety of complex quantum impurities and 
phase transitions, such as the interplay of Kondo effect 
and inter-impurity couplings [71-73], coupling to super¬ 
conducting or Dirac fermion baths [74-78] and the effect 
of multi-levels or multi-channels [79-82], see Refs. [24, 83] 
for a review. The fidelity susceptibility will provide a 
valuable tool to discover rich physical phenomena in such 
settings. In a broader context, since the Kondo effect in 
the quantum impurity models is often linked to the Fermi 
liquid behavior in the framework of DMFT, calculating 
the fidelity susceptibility Eq. (5) of the auxiliary quantum 
impurity problems may shed light on phase transitions of 
correlated materials. 
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